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Quantum dragon segments are nanodevices that have energy-independent total transmission of 
electrons. At the level of the single-band tight-binding model a nanodevice is viewed as a weighted 
undirected graph, with a vertex weight given by the on-site energy and the edge weight given by the 
tight-binding hopping parameter. A quantum dragon is a weighted undirected graph which when 
connected to idealized semi-infinite input and output leads, has the electron transmission probability 
T(E) = 1. The probability T(E) is obtained from the solution of the time-independent Schrodinger 
equation. A graph must have finely tuned tight-binding parameters in order to have T(E)= 1. This 
papers addresses which weighted graphs can be tuned, by adjusting a small fraction of the total 
weights, to be a quantum dragon. We prove that with proper tuning any nanodevice can be a 
quantum dragon. Three prescriptions are presented to tune a weighted graph into a nanodevice. 
The implications of the prescriptions for physical nanodevices is discussed. 

PACS numbers: 72.10.Bg, 73.23.Ad, 73.63.-b, 73.63.Fg 


INTRODUCTION 

Recently, the theoretical discovery of quantum dragons 
has been published [l|. This paper addresses the perva¬ 
siveness of quantum dragon nanodevices. 

All who are taught the basics of electricity learn about 
electrical resistance. Electrical resistance is the physical 
principle underlying the glow of light bulbs, the heat¬ 
ing of the element in your conventional oven and of your 
stovetop, the heat produced by your laptop computers 
and mobile devices, and the energy losses in transmission 
lines. A perfect electrical conductor would have zero elec¬ 
trical resistance, and consequently classically an electric 
current would not cause any heating. 

A recent discovery of a large class of nanodevices called 
quantum dragons [l| has potential technological impacts 
in nano-electronics. A quantum dragon is a nanode¬ 
vice that has zero electrical resistance, and therefore is 
a perfect conductor. The caveat is that in the quantum 
regime the electrons propagate coherently, and for co¬ 
herent propagation the classical concept of Ohm’s law is 
not valid |2|-[8[]. A quantum dragon has zero electrical 
resistance in a four-probe measurement [l|. In a two- 
probe measurement, the electrical resistance of a quan¬ 
tum dragon is quantitized. In particular, for a single open 
quantum channel the electrical resistance in a two-probe 
measurement is Gq 1 , with the quantum of conductance 
Go = 2 e 2 /h — 7.748 x 10 -5 S [3-8]. Here e is the charge 
of an electron and h is Planck’s constant. 

The seminal work of Landauer 0 showed that from 
the solution of time-independent Schrodinger equation 
0 the transmission probability as a function of energy, 
T(E ), was critical to obtain the electrical conductance 
via what is now called the Landauer formula. The Lan¬ 
dauer formula is for two-probe measurements. The elec¬ 


trical conductance is given by a convolution of the Fermi 
function for the electrons with T{E). For more compli¬ 
cated measurements, such as a four-probe measurement, 
the Buttiker-Landauer formula is the generalization of 
the Landauer formula that must be used. Many excellent 
books have been written using the Buttiker-Landauer for¬ 
mula [3^8J, including coherent transport through nano¬ 
structures 

A quantum dragon is a nanodevice that when con¬ 
nected to idealized leads (as in the example of Fig. [1]) 
has full transmission of electrons for all energies that 
can propagate through the leads [l[. In other words, 
in a quantum dragon electrons for any propagating en¬ 
ergy E have total transmission, the transmission prob¬ 
ability is T(E)= 1. Thus a quantum dragon nanodevice 
will behave as if the electrons underwent ballistic propa¬ 
gation. However, ballistic propagation requires no scat¬ 
tering of incident electrons, while there may be strong 
scattering of electrons in a quantum dragon. The best 
known experimental example of a quantum dragon is a 
single-walled carbon nanotube in the armchair configu¬ 
ration [IM3. For a quantum dragon T(E)= 1, and in a 
two-probe measurement the Landauer formula thus gives 
the single-channel electrical conductance Go- 

At the most basic level, the single-band tight-binding 
model for electrical transmission can be viewed as a 
weighted undirected graph [H, [23] connected to idealized 
semi-infinite leads. The tight-binding parameters include 
the on-site energy associated with a the vertex weight 
(atom weight or node weight) of the graph. Vertex j will 
have a weight 6j. The tight-binding hopping parameters 
are the weights given to the edges (bonds) of the graph. 
We consider graphs composed of £>1 slices, with only 
electron hopping allowed within a slice or between atoms 
in nearest-neighbor slices. We use tjj/ as the strength 
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of the intra-slice hopping parameters between vertices j 
and j'. Similarly, we use Sjj> as the strength of the inter¬ 
slice hopping parameters between vertices. Similarly, Wj 
or Uj are the strengths of the weights of hopping between 
the last (first) vertex of the input (output) lead and each 
vertex in the first (last) slice. The statements and proofs 
in this paper are in regard to weighted undirected graphs. 
The discussion section will bring the weighted undirected 
graph study back to examination of physical nanodevices. 

In this paper, we show the weighted undirected graph 
associated with any nanodevice can be a quantum 
dragon. The question to be addressed is how many pa¬ 
rameters within the tight-binding model must be tuned 
in order to obtain T(E)= 1. Three prescriptions are 
presented to make the weighted graphs into a quantum 
dragon. In the first two prescriptions, we consider only 
homogeneous weighted graphs made of i identical slices 
of m atoms, with only the simplest inter-slice connec¬ 
tions. In the third prescription the constraint of homo¬ 
geneity of slices is relaxed. The three prescriptions to 
make a weighted undirected graph a quantum dragon by 
tuning certain tight-binding parameters are listed below. 
In all cases we only consider connected graphs for the 
nanodevice, although some of these prescriptions may be 
modified to deal with graphs that are not connected. 

1. If all Ci and intra-slice hopping tjj/ are arbitrarily 
fixed, and only the simplest inter-slice hopping pa¬ 
rameters Sj jf are included, a weighted graph com¬ 
posed of identical slices can be a quantum dragon 
by tuning 2m+2 parameters: 

• The hopping parameters between the lead and 
nanodevice must each be a unique value; 

• A uniform electrical potential shift of all on¬ 
site energies ej] 

• All inter-slice hopping parameters Sjj/ must 
be scaled (tuned) by the same value; 

• The output lead connections to the last (£ th ) 
slice are identical to the input lead connections 
to the first slice. 

2. If all input lead to nanodevice hopping parameters, 
Wj, are in arbitrarily fixed ratios and non-zero, only 
the simplest inter-slice hopping terms are included, 
and all €j and intra-slice hopping tjj* are fixed, 
requires tuning 2ra+2 parameters: 

• Each ej must be shifted (tuned) by a unique 
value of the electrical potential at site j by Vj ; 

• The strength of all lead-device hopping pa¬ 
rameters must be normalized (tuned) in a 
unique fashion; 

• All inter-slice hopping parameters must 
be scaled (tuned) by the same value; 


• The output lead connections, Uj, to the last 
(^ th ) slice are identical to the input lead con¬ 
nections, Wj , to the first slice. 

3. If all ej and intra-slice hopping tjjf are arbitrarily 
fixed in each slice, but the slices of the weighted 
graph are not homogeneous: 

• The hopping parameters between the leads 
and the nanodevice must each be a unique 
(tuned) value, which may be different for the 
connections to the input lead ( Wj ) and output 
lead {uj)\ 

• A uniform electrical potential shift (tuning) 
of all on-site energies ej in each slice, but the 
tuning required may be different for each slice; 

• A unique (tuned) value for each of the inter¬ 
slice hopping terms Sjjt, which may be dif¬ 
ferent between each pair of atoms on adjacent 
slices. 

The vast majority of calculations use the Green’s func¬ 
tion method to calculate the electrical conductance [3|-[8j . 
The Green’s function implicitly solves the Schrodinger 
equation, and through the Landauer equation !2] gives 
T(E) and the electrical conductivity. In this paper, 
we use a different standard technique, called the ma¬ 
trix method H3 , which solves the Schrodinger equation 
using an particular ansatz. From the matrix method 
solution of the Schrodinger equation, one can calculate 
T(F), and hence the electrical conductivity from the 
Landauer equation [2]. There are several papers that 
use the matrix method to calculate T ( E ) 0,013. The 
physical quantity, T(E ), is independent of which calcu¬ 
lation method is used since they both satisfy the time- 
independent Schrodinger equation for the semi-infinite 
leads connected to the nanodevice. 


TIGHT-BINDING MODEL 

A simple example of a weighted undirected graph asso¬ 
ciated with a physically realizable nanodevice that may 
be made into a quantum dragon is shown in Fig. [U A 
simple-cubic (SC) lattice structure is presented in Fig.[TJ 
The metal polonium (Po) has a simple cubic lattice struc¬ 
ture, and hence Fig. [l] can be viewed as a nano-crystal of 
Po. Fig.QJa) shows a single slice of the nanodevice, made 
up of all atomic sites in a square lattice that fit into a 
given ellipse. The center of the ellipse is randomly given 
within a square lattice unit cell, leading to a fixed arbi¬ 
trary non-isotropic atomic arrangement. The first slice 
of the nanodevice is connected to the input lead (wj), 
tuned as shown in Fig. QJb). The last slice of the nan¬ 
odevice is connected to the output lead in an identical 
fashion, namely Uj=Wj for all j. An example of a twenty 
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With these units, only electrons with energies —2 < E < 
2 propagate in the leads mini- 

For the tight-binding model the hopping parame¬ 
ters are obtained from the kinetic energy terms of the 
Schrodinger equation [3], and therefore must be non¬ 
positive. That is the reason for the negative sign in the 
hopping parameter —1 within the leads. We will take 
all hopping parameters ( tjj /, Sjjt , Wj , and Uj) to be 
positive, i.e. the hopping strengths are positive. The re¬ 
quired negative signs are put in explicitly, so tjj />0 and 

s j,j' >0. 

TRANSMISSION VIA THE MATRIX METHOD 



FIG. 1. (Color online.) See the Appendix for full information 
on the figure construction, (a) An example of a slice used 
to make a quantum dragon, here composed of atoms from 
a square lattice within an ellipse (green curve), with nearest 
and next-nearest neighbor bonds (blue cylinders, which rep¬ 
resent the intra-slice hopping terms tjj/). There are m =76 
atoms and 236 bonds in the slice, (b) The hopping terms Wj 
(orange cylinders) between a lead atom (white sphere) and 
the first slice of the nanodevice are shown. As discussed in 
the text, the lead-device hopping terms must have particular 
values (denoted by the radius of the orange cylinders) for the 
nanodevice to be a quantum dragon. The lead atom is places 
above the first slice, and located at the center of mass (CM) of 
the connections Wj. (c) The complete nanodevice composed 
of ^=20 identical slices, connected to the two semi-infinite 
leads (only two atoms in each lead are shown). 


slice (£=20) nanodevice connected to the leads is shown 
in Fig. (He), one which may be tuned to be a quantum 
dragon. 

We assume that the on-site energy of the semi-infinite 
leads is zero, thereby setting our zero of energy. We 
take the hopping parameters between the lead atoms as 
— 1, thereby setting our energy scale as the strength of 
the lead-lead hopping term. We take the lattice spacing 
within the leads to be unity, setting the unit for length. 


In order the calculate the transmission of an incom¬ 
ing electron as a function of energy, T(E ), the time- 
independent Schrodinger equation needs to be solved. 
The matrix method for solving the tight-binding model 
is used here [lg]. This method has been used and pub¬ 
lished a sufficient number of times that the derivation of 
the equations [1] is not required. Therefore, only the rel¬ 
evant equations are given in order to set the notation [l| 
and give the reader the basics of the matrix method. 

The nanodevice is composed of £ slices. For the first 
two quantum dragon prescriptions, the slices are identi¬ 
cal, each with m atoms. The equation to solve for the 
Schrodinger equation of the device (with £m vertices) and 
the semi-infinite leads is infinite. After using an ansatz 
for the leads, the final matrix equation has a linear di¬ 
mension £m +2 0,0. We assume that only atoms in 
nearest-neighbor slices interact. Fig. [T] shows an exam¬ 
ple with m =76 and ^=20. For £=4 slices, the matrix 
equation to analyze has the form 
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of coupling the finite graph associated with the nanode¬ 
vice to the semi-infinite leads. The zero matrices, 0 , are 
mxm, and the zero vectors 0 have m elements. The 
inter-slice hopping matrix B is mxm and has elements 
of the inter-slice hopping parameters. For the simplest 
inter-slice connections, as shown in Fig. [U^c), the inter¬ 
slice hopping matrix is given by B = — sq! where I is the 
mxm identity matrix. The vector w has as its elements 
the hopping parameters between the last atom of the in¬ 
put lead and the atoms in the first slice. The vector u 
has as its elements the hopping parameters between the 
first atom of the output lead and the atoms in the last 
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(£ th ) slice. The mxm matrix F = A — El, with E the 
energy of the incoming electron. The mxm matrix A 
has as its j th diagonal element the on-site energy ej of 
the j th atom (vertex). The j,j' off-diagonal element of 
A is the intra-slice hopping term —tjj/. The matrix A 
is symmetric, since the Schrodinger equation involves a 
Hamiltonian, and we here restrict ourselves to real values 
for the on-site energies and all the hopping parameters. 
In Fig. [IJa) the intra-slice hopping terms are shown as 
cylinders, with the radius of the cylinder representing the 
strength of the intra-slice hopping term tjjr . 

The matrix equation to solve for the transmission is 
given by (again written for £=4 slices) 


l + r \ 
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with the definition ’B.(E) = — 2 iXm{^{E)) [l, 16]. The 
wavefunction of slice j is given by ijjj . For a given energy 
of the incoming electron, the inverse of the matrix of N 4 
in Eq. © is calculated. This enables one to obtain the 
wavefunctions ipj, as well as r and tr- Note that these 
all depend on the energy E of the incoming electron. 
The quantity tr should not be confused with the hopping 
parameters in N^; unfortunately both are conventionally 
denoted by a lower case £, so the subscript T denotes that 
tr is for the transmission, not a hopping parameter. The 
transmission probability of the electron is given by 

T(E) = \tr(E) | 2 . (3) 

The reflection probability for an electron is given by 
7Z=\r\ 2 . Every electron is either reflected or transmit¬ 
ted, so 7£+T=l. 


with the mxm matrices X unitary [i.e. X^X=I). Note 
that (written for £=3) 
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( 5 ) 


Form the matrix M< = X^NgXj. In the standard fash¬ 
ion, by matrix multiplication the mxm diagonal matrix 
blocks in Mf are given by XFXT The non-zero off- 
diagonal blocks of M/ : are given by XBX^ and XT. or 
by the complex congugates. For example, for £=3 


M 3 


/ £(E) #Xt 0+ 

XT XFXt XBXt 
0 XB+Xt XFXt 
0 0 XBtXt 

\0 0 + Qt 


0 f 


0 

XBXt 

XFXt 

fftxt 


° \ 
0 

0 


Xu 

m) 


■ (6) 


Multiplying Eq. <[2|) on the left by X(, inserting XjX^==I 
between and the vector in Eq. © that contains the 
wavefunctions, and using Eq. ([5|) gives (written for £=4) 
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We have complete freedom to choose the unitary trans¬ 
formation matrices X. Assume we can find a X that 
satisfies the four mapping equations 0 ] 


XAX ' = ( 5 A ) ■ (8) 

xbx ' = (T t ) ■ 


PROJECTION MAPPING TECHNIQUE / ~ \ 

X ^ = (T)’ (10) 

In order to more easily calculate tx(E) in Eq. © we 
introduce a transformation matrix X, written for (=3, and 
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( 11 ) 


(4) The (m— l)x (m—1) matrices A and B are not important, 
since they will not be connected by any path to either the 
input or output leads [ 1 ]. 
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Introduce the (£+ 2 )x(£+ 2 ) matrix Mg, written for 
£=4, as 
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with k = e — E. This is the matrix formed from only 
the £ transformed-sites in Eq. 0 that are connected to 
the leads, after the mapping equations Eq. (j 8 ]) through 
Eq. ITT]) are used. 

The probability of transmission of the electron of en¬ 
ergy E, T=\tr\ 2 are calculated from tr found from either 
the equation (written for ^=4, with the displayed vectors 
of length 4ra+2) 
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(13) 


or from the equation (written for ^=4, with the displayed 
vectors of length 6=4+2) 
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where the <j)j are the first element of the vector X 0 j. 
The mapping equations of Eq. (J 8 ]) through ([TT]) thus re¬ 
duce substantially the size of the matrix that one must 
find the inverse of in order to calculate tr- The_ma- 
trix is (£m+2) x (£m+2) while the matrix is 
(£+ 2 )x(£+ 2 ). Note that provided the mapping equa¬ 
tions [Eq. (j 8 j) through ([TT]) ] hold, no approximation is 
made by going from Eq. (fT3]) to Eq. m , as they both 
give identical transmissions T(E)=\tr \ 2 . 


QUANTUM DRAGONS FROM MAPPING 

Recall that our lead sites have on-site energies set to 
zero and a hopping of strength unity. The nanodevice 
will be a quantum dragon, i.e. will have 7 ~(E) = |£t| 2 =1, 
if e=0 and s w =s u =§b= 1. The reason for the complete 
transmission of electrons of all energies is that these are 
the values that the matrix in Eq. m would have if 
one used the matrix method to calculate the transmis¬ 
sion through a homogeneous infinite wire, but selected a 


string of i lead sites to be a nanodevice. Since the wire is 
homogeneous, there is no scattering by the nanodevice, 
and the electrons for any energy that propagate through 
the lead would be completely transmitted. Note, how¬ 
ever, that the original slices (as in Fig. [T]) may be very 
inhomogeneous, which may lead to strong scattering of 
the electrons. 

There is complete freedom in terms of the transforma¬ 
tion matrices X that are used. The only requirement 
to go from Eq. CG2> to Eq. ([T4]) is that the four map¬ 
ping equations, Eq. (J 8 ]) through ([TT]) . are satisfied. For 
the first two quantum dragon prescriptions, only nanode¬ 
vices with the simplest inter-slice coupling between iden¬ 
tical slices (as in Fig.[]Jc)) will be analyzed. This means 
the inter-slice hopping matrix is B=— 85 I. Hence Eq. (J9]) 
is satisfied since X is unitary. Thus one has for these 
two prescriptions 55 = 85 . In order to have the nanode¬ 
vice have the possibility of being a quantum dragon we 
tune the parameter 85 = 1 . Therefore, for such a quantum 
dragon nanodevice we need only to find a transformation 
matrix X that satisfies Eqs. ([ 8 ]), ITOl) and ([TT]) . and that 
have e =0 and 8 w =s u =l. 

The third quantum dragon prescription will require the 
four mapping equations, Eqs. ([ 8 ]) through ([TT]) to be sat¬ 
isfied, and to have for each slice k the mapped on-site 
energy e/c= 0 , between each nearest-neighbor pair of slices 
5/e,/c+i=l, between the input lead and the first slice s w = 1 , 
and between the output lead and the last (£ tYl ) slice s u = 1 . 


QUANTUM DRAGONS: PRESCRIPTION 1 

We consider a weighted undirected graph made from 
i identical slices each with m atoms. As in Fig. Q] the 
weighted graph may be associated with the tight-binding 
model on a physical nanodevice. We assume the intra¬ 
slice matrix elements A are all fixed to arbitrary values 
(both the on-site energies ej and intra-slice hopping terms 
tjj'). We assume that the slice atoms are strongly con¬ 
nected, in the sense of graph theory [Hi [HI- In other 
words, within a slice every atom can be visited starting 
from any other atom by a series of hops using only the 
non-zero intra-slice hopping terms tjj /. We assume that 
only the simplest inter-slice connections are present, so 
B = — 80 1. We assume we are free to tune all the lead- 
slice connection strengths in w and u, to add a constant 
electric potential V to the on-site energy of every slice 
(which is the same for every atom in the nanodevice), and 
to tune the strength of the inter-slice hopping strengths 
8 0 . 

Let e max be the maximum of zero or the largest positive 
diagonal element of A, i.e. 

£max = max{0, ei, e 2 , • • • , e m } (15) 

with the €j the on-site energies, which are the diagonal 
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elements of the matrix A. Introduce the matrix 


A ^rnaxl A. 


(16) 


Then the matrix A is non-negative 24. 25j|. which is writ¬ 
ten as A > 0. In other words, every element of A is 
positive or zero. Let w\ be the normalized eigenvector of 
A associated with the largest eigenvalue Ax of A, i.e. 


Aw i = Xiwi 


(17) 


Since A has an associated strongly connected graph 
and is non-negative, by a well-known extension of the 
Perron-Frobenius theorem snni the vector w\ is unique 
and can be written with all non-negative values (Fq>0). 
Then w\ is also an eigenvector of A with eigenvalue Ai, 
i.e. 


Awi — Ai w\ . (18) 

Remember that the eigenvalues are functions of all the 
elements of A, including the largest diagonal element of 
A, i.e. 


Al — Ax (c max ) — C m ax Ax (Cmax) • (19) 

Note that A and A are symmetric, so that they can be 
associated with an undirected graph, as opposed to being 
associated with a directed graph as would be required for 
general non-negative matrices. Since A is symmetric, w\ 
is both a right and left eigenvector of A. 

We tune the lead-device hopping terms, which must all 
be non-positive, to be 

w = —s w w\ and u = —s w w\ (20) 

with some overall strength s w that we will determine be¬ 
low. The strengths are the same for the connections to 
the input and to the output leads. We have had to form 
the matrix A in order to use the Perron-Frobenius the¬ 
orem to show that one can obtain vectors w and u with 
all non-positive values. This is the physical constraint 
imposed by the hopping being the negative of the ki¬ 
netic energy portion of the time-independent Schrodinger 
equation. 

Choose the transformation matrix to be (written for 
771=5) 


= ( W\ W 2 Ws W 4 W$ ) (21) 

with orthonormal vectors WjWj/=5jj/. The vector w\ is 
an eigenvector of A from Eq. m • The other vectors 
Wj for j=2, 3, • • • ,777 need not be eigenvectors of A, only 
orthonormal to each other and to w\. With this choice 
of X and w, the mapping equation Eq. m is satisfied 
with s w =s w . Similarly, the mapping equation Eq. m is 
satisfied with s u =s w . The connections to the input and 
output leads are the same since all slices are identical. 


We tune one parameter by applying a constant electri¬ 
cal potential — Ehift to every atom in the nanodevice, so 
the matrix A is shifted to the matrix A—EhiftI- Note 
the electric potential shift is the same for all atoms in the 
nanodevice. Then 


(A - EshiftI) wi = (Ax - Ehift) wi • (22) 

With our choice of X, the mapping equation Eq. m is 
satisfied with 


e — Ai(e max ) Eshift • (23) 

All four mapping equations have thus been satisfied. 
For a quantum dragon using this prescription, we need 
only choose 


w = —w\ 

u = —W\ 

and 

E s hift = Ax 


(24) 


so that s w =s u =1 and e=0. We have tuned the m lead- 
slice interactions in w = —tFx, the m lead-slice interac¬ 
tions in u = — tFx, the constant shift electrical potential 
Ehift? an d the inter-slice hopping strength s^m 1. In all, 
we have tuned 2t?7+2 tight-binding parameters. With the 
simplest inter-slice connections B=—sqI? the total num¬ 
ber of possible parameters in the tight-binding model is 
equal to 


Axotal 


parameters 


+ 3m + 1 

m 2 +5m+2 

2 


(25) 


The number of intra-slice hopping terms tjjf is 777(777 — 
l)/2, the number of different on-site energies ej is 777 , the 
number of lead-slice hopping terms is m for the input 
lead and m for the output lead, and there is the inter¬ 
slice hopping strength so, giving the result in Eq. ([25]) . 

This prescription does not require either that there be 
any symmetry within a slice or that the underlying graph 
for a slice is planar. Most importantly, we have not had to 
tune the connections or the intra-slice hopping strengths 
in a slice. The slice may be semi-regular, as in Fig. [1] 
but that is not necessary. One could have each slice, 
for example, be a portion of a two-dimensional quasi¬ 
crystal. One can also have the slices atomic arrangement 
be completely random, as in an amorphous material. One 
example of such a complicated, amorphous nanodevice 
that can be a quantum dragon is shown in Fig. [2] In 
this prescription, quantum dragons exist everywhere on 
a 2 777 + 2 dimensional ‘surface’ in the ( 777 2 + 5tt 7 + 2)/2 
dimensional parameter space. 


QUANTUM DRAGONS: PRESCRIPTION 2 

In the second prescription, we again assume that all 
slices are identical, and that only the simplest inter-slice 
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FIG. 2. (Color online.) See the Appendix for full informa¬ 
tion on the figure construction, (a) A slice of a nanodevice 
made from random arrangements of atoms inside of an ellipse 
(green). Here there are m —50 atoms in the slice, (b) The 
lead-device hopping terms Wj required to make the nanode¬ 
vice comprised of slices as in (a) into a quantum dragon. The 
radii of the lead-slice (orange) cylinders are proportional to 
the strength of the required hopping terms. The positioning 
of the lead atom is at the CM of the Wj bonds, but displaced 
above the slice along the direction of electron propagation 
in the leads, (c) The complete nanodevice formed from £=8 
identical slices, connected to the lead atoms (white spheres) 
via the lead-slice hopping terms gives a quantum dragon. 


hopping is present (B=—sol)- We further assume we 
have fixed arbitrary lead-device connections wo, and any 
fixed arbitrary intra-slice matrix Ao- The subscript 0 
stands for the original given problem. We assume every 
element of wo is negative, (— FJo>0). We are required to 
keep fixed the ratios of the elements of but can tune 
the overall normalization. Therefore the final lead-site 
hopping connections are given by 

w = s w wq (26) 

where we can only tune s w . Similarly, the connection to 
the output lead is given by 

u = s u w 0 (27) 

where we can only tune s u . Choose a coordinate sys¬ 
tem so the slice-to-slice direction of electron propagation 
is along the z-axis. We assume that the graph associ¬ 
ated with Ao is strongly connected, even though Ao may 
have many elements equal to zero. We further assume we 
are only allowed to tune the intra-slice matrix by adding 
an electrical potential V(x,y) at the location of every 


atom. The electric field is given in the standard fash¬ 
ion, E(x,y) = —WV(x,y). The electric potential must 
be continuous, but otherwise its values at points between 
atoms does not enter the tight binding model. The elec¬ 
tric field at the device location is the same for all slices. If 
atom j has the coordinates (xj,yj), the j th diagonal ele¬ 
ment of Ao is thus changed by tuning the term V{xj , yj ). 
Introduce the mxm diagonal matrix Dy with the j th 
diagonal element equal to the tuned electrical potential 
values, V(xj,yj). The new intra-slice matrix with this 
added electrical potential is 


A = Ao + Dy . (28) 

Changing the electrical potential has not changed the 
intra-slice hopping strengths (tjjr ) or the graph connec¬ 
tivity of a slice. 

We choose the electric potential so that —A > 0, i.e. 
so all elements of A are non-positive. We want to choose 
a Aq and the m diagonal elements of Dy so that 


(-Dy - A 0 ) wo = A 0 FJ 0 (29) 


and — (Dy+ Ao)>0. Introduce the vector e with all 
elements unity. Since the given wq has no zero elements, 
there is a diagonal matrix D e such that D e FJo «= e and 
the inverse (diagonal) matrix D” 1 exists. Also introduce 
the square matrix J with all elements unity. One has 
J=eet, and for the J matrix mxm one has —Je=e. The 
diagonal matrix Dy we need to satiisfy Eq. ([29]) is 

D v = —A 0 I - IO (DeAoD; 1 J) (30) 

where O is the Hadamard (element-by-element) matrix 
product. The proof is that 


DyFJo 


[—AoI-lO (DeAoD-iJ)] «?o 

-\ 0 W 0 - D e (10 (AoDyee 1 ')) w 0 
-X 0 w 0 - D e (I © (AotooO 
-Aotoo - (I © (A 0 MJ 0 e 1 ')) D e w 0 
-A 0 «J 0 - (I © (AouJoO e 

— XqWo — AqWo 


as required in order to satisfy Eq. ([29]) . 

Introduce the magnitude of wq to be |hJo| = xjW qWq. 
Now choose the transformation matrix X to have the 
form as in Eq. m , with w\ - —wo/ |reo|, again with or¬ 
thonormal vectors w\wj = for = 1,2, ■ ,m. 

With this choice for the transformation matrix X the 
mapping equation Eq. m is satisfied with s w = \wo \ s w . 
Similarly, the mapping equation Eq. (HU is satisfied with 
s u = |rco| s u . The lead-device connections will allow a 
quantum dragon if one chooses s w =s u =l, which means 
we need to choose s w =s u = 1/ |reo|. 

We have the simplest possible inter-slice hopping 
terms, having chosen B = — sqI- As in prescription 1, 
since XX^=I the mapping equation Eq. m is satisfied. 
The inter-slice terms allow a quantum dragon if sq=1. 
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With the tuned on-site energies from Dy of Eq. ([30]) . 
the mapping equation Eq. (|8j) is satisfied with e=—Ao- 
We are free to tune Ao, since this would involve an equal 
shift of the electric potential on every site of the nan¬ 
odevice. In order to have the intra-slice matrix allow a 
quantum dragon requires e= 0, and therefore we tune to 
make Ao=0. From Eq. m the required electric potential 
that must exist at each atomic site is 

Dy = - I © (DeAoDj 1 J) . (32) 

Therefore, we have found the solution to all four map¬ 
ping equations in order to have a quantum dragon, i.e. to 
have T(E)= 1, for prescription 2. We have tuned 2m+2 
tight binding parameters in order to have a quantum 
dragon, from the (m 2 + 5 m + 2)/2 total parameters [see 
Eq ([25]) ] . An example of the tuning of each slice in pre¬ 
scription 2 is illustrated in Fig. [3] 

In both prescriptions 1 and 2 the inter-slice mapping 
equation, Eq. m, has been satisfied with the inter-slice 
matrix B = —sqI- However, Eq. could also be satis¬ 
fied by having an inter-slice matrix of the form 

B = —(33) 

W{ W\ 

In order to have a quantum dragon from the mapping, 
Eq. ([33]) must be satisfied and we must tune the hopping 
strengths to satisfy So+s ww =l. The third prescription is 
a generalization of this tuning of the inter-slice hopping 
terms to the case where every slice may be different. 


QUANTUM DRAGONS: PRESCRIPTION 3; 

INHOMOGENEOUS SLICES 

We make similar assumptions about the intra-slice in¬ 
teractions as in prescription 1. However, we now let every 
slice possibly be different. We assume the intra-slice ma¬ 
trices A/e for slices 1 <k<£ are given fixed values, and all 
slices are assumed to be strongly connected. We assume 
we are again allowed to tune a constant electric potential 
V k by a constant value for every atom of slice k. We as¬ 
sume we are allowed to tune the lead-device connection 
vector w for the input lead, and the lead-device connec¬ 
tion vector u for the output lead. We further assume 
that we are allowed to tune every element of the inter¬ 
slice hopping matrix B/^+i for k= 1, 2, • • • (£ — 1). 

For this non-identical slice nanodevice, the matrix 
equation to analyze has the same form as Eq. with 
two differences. The diagonal matrices (intra-slice terms) 
F k =A k —El may all be different. The matrices F k are 
m k xm k , where slice k has rrik atoms (vertices). The 
inter-slice terms, B/^/c+i may all be different, and are in 
general not square matrices, being of size mjtxmqi. In 
general, the device-output-lead interaction u will usually 
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FIG. 3. (Color online.) (a) A slice of a nanodevice made from 
square lattice arrangements of atoms inside of a generalized 
ellipse and a circle. There are ra=60 atoms in the slice, (b) 
A lead-slice interaction connection, chosen without regard to 
the arrangement of atoms or bonds in the slice. The (black) 
sphere is the last atom of the input lead and the (black) bonds 
are the hopping parameters (given by the vector wo) between 
this lead atom and the atoms in the first slice. The strength 
of the lead-slice interaction is proportional to the radius of 
the (black) cylinders. The location at which the lead atom 
would have to be placed in order to be at the CM of the 
wo is shown as the yellow sphere, (c) The required electrical 
potential at each atomic site to make the device a quantum 
dragon. The height of each (cyan) cuboid is proportional to 
the value of V at the atomic site to make T(E)= 1. For visual 
reasons, all potentials are shifted by a constant to make them 
all non-negative. 


be tuned to be different from the device-input-lead inter¬ 
action w. 

For our non-identical slice case, to perform the map¬ 
ping the transformation matrix X of Eq. will have 
different unitary transformation matrices X/~ of size 
m k xm k along the diagonal. The four mapping equations 
of Eq. (|8j) through Eq. (fill) now become a set of 2^+1 
mapping equations. The £ intra-slice mapping equations 


















are 


XfcAfcXj. = 


h O f \ 

0 A f. ) 


The £—1 inter-slice mapping equations are 


X fc B 


k*5k,k+l 


yT _ 
■ A fc+1 


—h,k 0 f A 

o B k ) 


of Afc. Since every A k has an associated strongly con¬ 
nected graph and is non-negative, by the extension of 

( 34 ) the Perron-Frobenius theorem 24|, 25] every vector w\ :k 
is unique and can be written with all non-negative values 
(uJi,*>0). Hence, w\,k is also an eigenvector of A k with 
eigenvalue Aq^. 

We tune the lead-device hopping terms for the incom- 

(35) ing lead, which must all be non-positive, to be 


The lead-device hopping terms for the input lead must 
satisfy the mapping equation 


Xi w = 



(36) 


w = -s w wi, 1 (41) 

with some overall strength s w . Similarly, we tune the 
lead-device hopping terms for the outgoing lead, which 
must be non-positive, to be 


while the lead-slice hopping terms for the output lead 
must satisfy the mapping equation 

=(T)' (37) 


The (ra*; —l)x(rafc — 1 ) matrices A k and the 
(rrik — 1 ) x (m k +i — 1 ) matrices are not im¬ 

portant, since they will not be connected by any path to 
either the input or output leads [l|. 

As for Eq. (HD, for the non-identical slice case we in¬ 
troduce the (^+2)x(^+2) matrix M^, written for £=4, 
as 


M 4 


(S(E) -s w 

0 -s b , l 
0 0 
0 0 
\ 0 0 


0 0 

-Sb, i 0 

^2 —Sb, 2 

— $b,2 

0 —Sb, 3 

0 0 


0 0 \ 

0 0 

0 0 

—h, 3 0 

/^4 S u 

—s u £{E) ) 


(38) 


with kk = e k —E. This is the matrix formed from only 
the £ transformed-sites after the 2^+1 mapping equations 
Eq. m through Eq. ([37]) are used. The probability of 
transmission of the electron T = |£t | 2 is calculated from 
the tT quantities found from the solution of the matrix 
Eq. (1T41) . 

To satisfy the intra-slice mapping equations in 
Eq. ([34]) . follow the same procedure based on the same 
arguments as in Eq. ([T5]) to Eq. ([241) . In particular, intro¬ 
duce the maximum of zero or the largest positive diagonal 
element of A&, i.e. 


e max,/c — max{0, * * * 5 e mfc,/c} (39) 

with the €j : k the on-site energy of the j th atom in slice 
k. Introduce the matrix 

Afc — £max,fcl Afc , (40) 

and hence A k > 0. Let wi, k be the normalized eigen¬ 
vector of Afc associated with the largest eigenvalue Ai^ 


u = -s u wi,£ (42) 

with some overall strength s u . We have thus far had to 
form the matrices Ai and A^ in order to use the Perron- 
Frobenius theorem to show the two lead-device vectors w 
and u are non-positive, as required for a physical tight- 
binding model. 

Choose the transformation matrices to be (written for 

m k =6 ) 

X£ = ( Wi ik W 2 ,k ™3,k W4,k W6,k ) (43) 

with orthonormal vectors k Wj f ,k=5j,j' • The vector 
w± jk is an eigenvector of A k . The other m k — 1 vectors 
w jik for j= 2,3, • • • , m k need not be eigenvectors of the 
A*;, only ortho normal to each other and to wi jk . With 
this choice of X& and the tuned values for the lead-device 
vectors w and u, the mapping equations ([36] ) and ([37] ) are 
satisfied with s w =s w and s u =s w . 

We apply a constant electrical potential — V^hift ,k to 
every atom in the slice k of the nanodevice, so every 
matrix A k is shifted to the matrix A k — V^hift,fcl- Then 

(A k Fshift,/cl) W lt k = (Al,fc Fshift,/c) ^1 ,k • (44) 

With our choice of X/~, the mapping equation Eq. dm is 
satisfied with 


e k — Ai 5 / C (e max ,/ C ) Fshift,/c • (45) 

In order to satisfy the mapping equations in Eq. ([35]) 
we tune the inter-slice hopping matrices, which are 
rafcxrafc+i, to be 

Bfe, fc +i = -s b)k wi )k w\ k+1 (46) 

with arbitrary inter-slice hopping strengths s b , k . We have 
had to form the matrices A k in order to use the Perron- 
Frobenius theorem to show that all elements of ’B kjk +% 
are non-positive. 

All 2^+1 mapping equations for the non-identical slice 
case have thus been satisfied. For a quantum dragon 
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using this prescription, we need only further tune the 
values so 


so that s w = 1, 
l<k<L 


w 

= -^1,1 , 


u 

= -wi,e 


bshift,/c 

= Ay/c 

and 

(47) 

^k,k- (-1 

= -Wi, k w\ k+1 


11 

Sb,k= 1, and 4=0 

for all slices 


Each intra-slice matrix A& is symmetric, and hence has 
m k (m k + l)/2 tight-binding parameters. We have had 
to tune the diagonal m k of these by the same amount 
V^hiftso each slice has the one parameter T^hift ,k that 
needs to be tuned. There is only a unique choice for all 
inter-slice hopping terms ~Bk,k+i and the two lead-slice 
hopping terms w and u. The free tight-binding parame¬ 
ters are only from the intra-slice matrices A&, for a total 
number of arbitrary parameters 


A"free 


E 

E 


t 

k=l 

l 

k= 1 


m k (m k + 1) _ 
2 

(m fc +2)(m fc —1) 


2 


(48) 


The total number of tight-binding parameters is 


AT i TYlk ( mk + 1) , i TACi\ 

ATot = mi + ) J ---h> v m k m k +i + mi. (49) 

k=l k= 1 

An example of the formation of a quantum dragon using 
prescription 3 is shown in Fig.|U Prescription 3 does not 
require any translational symmetry, either along the di¬ 
rection of electron motion in the leads or perpendicular 
to this axis. Therefore, the nanodevice can be considered 
to be completely disordered. Nevertheless, all incoming 
electrons will be completely transmitted through the nan¬ 
odevice, i.e. it is a quantum dragon having T(E)= 1. 
The quantum dragons exist only on a low-dimensional 
‘surface’ [dimension given by Eq. ([48]) ] in the entire tight- 
binding parameter space [dimension given by Eq. ([49]) ]. 


DISCUSSION AND CONCLUSIONS 

We have shown that quantum dragons are ubiquitous. 
They exist for any fixed atomic bonding arrangement! 
We have chosen to concentrate on nanodevices with £ 
slices, but the prescriptions also work when there is only 
i—\ slice. For £=1, any atomic bonding arrangement 
between the atoms is possible. The only question in 
all prescriptions is how many tight-binding parameters 
need to be tuned, and to what values these parameters 
must be tuned. The prescribed tight-binding parameters 
must be tuned in order to satisfy the mapping equations 
Eq. ([8]) through Eq. ([IT]) . To allow the nanodevice to 
be a quantum dragon requires further tuning to specific 



r 



FIG. 4. (Color online.) A nanodevice which can be a quan¬ 
tum dragon. The device has a total of 44 atoms in £—A slices. 
The widths of the cylinders are proportional to the size of the 
hopping parameters. See the appendix for a full description, 
(a) Showing only the device atoms, lead atoms (two in each 
lead) and bonds, lead-device interactions, and intra-slice in¬ 
teractions. The lead atoms are positioned at the CM locations 
given by Wj and Uj for the first and last slice, respectively, (b) 
The figure in (a), with the inter-slice tuned hopping interac¬ 
tions (orange cylinders) in prescription 3 to make the device 
a quantum dragon. 


tight-binding values. With careful tuning, electrons of all 
energies will have complete transmission, T(A)=1, and 
the device will hence be a quantum dragon. 

Three prescriptions, all allowing for quantum dragons 
from inhomogeneous nanodevices, have been presented 
in detail. The first two prescriptions have the simplest 
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inter-slice hopping terms, and have every slice identical. 
The first prescription allows arbitrary fixed values for 
the positions of all atoms, of all electrical potentials (up 
to a constant shift), and of all the intra-slice hopping 
strengths between atoms. In the first prescription the 
lead-device hopping strengths must be tuned in a pre¬ 
scribed manner in order to obtain a quantum dragon. 
The second prescription is related to the first, but the 
lead-device connections are arbitrary values (but identi¬ 
cal connections for the input lead and output lead), while 
the electrical potentials on each atom must be tuned in a 
prescribed manner in order to obtain a quantum dragon. 
The third prescription has slices which may all be dif¬ 
ferent, with atomic bonding strengths and electrical po¬ 
tentials (up to a slice-dependent constant term) fixed ar¬ 
bitrarily, while the lead-device and inter-slice hopping 
terms must be tuned in a prescribed manner in order to 
obtain a quantum dragon. 

For all prescriptions the number of arbitrary parame¬ 
ters is much larger than the number of parameters which 
must be tuned in a particular fashion. For the first two 
prescription, the ratio of the number of tight-binding pa¬ 
rameters is 

tuned _ 4 (m + 1) , . 

Total m 2 + 5 m + 2 

where m is the number of atoms in every slice of the nan¬ 
odevice. For the third prescription the number of atoms 
rrik in each of the 1 <k<£ slices can be different. However, 
if all rrik=m (but the intra-slice bonds and electrical po¬ 
tential may be different for every slice) the ratio of the 
number of tight-binding parameters is 

tuned £m 2 + £m — 2t 

Total 3 £m 2 + £m — 2 m 2 + 4m ’ 

Thus in all cases, quantum dragons exist only on a 
low-dimensional ‘surface’ of the high-dimensional tight- 
binding parameter space. An analogy might be useful to 
understand the relationship between the complete tight- 
binding parameter space, the parameter space of the 
mapping equations, and the parameter space of quan¬ 
tum dragons. Consider a room, so the space has dimen¬ 
sion D= 3, which can be viewed as the complete param¬ 
eter space for this analogy. A thin sheet of paper in the 
room, maybe folded or crumpled, represents the param¬ 
eter space where the mapping equations hold, here D= 2. 
A curve drawn on the sheet of paper represents the pa¬ 
rameter space where quantum dragons exist, here D— 1. 
Clearly, a blind Monte Carlo search of the D =3 space 
would have zero probability of locating a point exactly 
on the D =2 surface, much less on the D =1 curve. 

The natural question is which, if any, prescription 
would yield a nanodevice and lead-slice connections that 
can be realized reasonably easily experimentally. The 
answer is that all three prescriptions have experimental 
difficulties. The first prescription can be used to have 


a single-slice (£= 1 ) nanodevice, and a quantum dragon 
can always be found with the correct lead-device connec¬ 
tions. However, to make the lead-slice connections would 
be physically impossible at the nanoscale, particularly 
for a non-planar arrangement of atoms in the slice. The 
second prescription for completely arbitrary lead-device 
connections would require electric fields on the order of 
10 9 V/m precisely tuned at the nanoscale level. Although 
not physically impossible, such high electric fields tuned 
to the nanoscale level would set an extremely high exper¬ 
imental bar. The third prescription requires inter-slice 
connections that seem impossible, even for the case of 
Fig. 0] with only about ten atoms per slice. 

A more realistic method of experimentally construct¬ 
ing an inhomogeneous quantum dragon might take an 
approach that is a combination of the three prescriptions. 
For example, one could require that the lead-device in¬ 
teractions should be monotonically dependent on the dis¬ 
tance from the lead atom, and the electric field required 
to change by no more than a few percent on the nanoscale 
level. These types of smoothness constraints are not nec¬ 
essary mathematically, but will be critical to synthesizing 
an experimental realization of inhomogeneous quantum 
dragons. With such smoothness constraints, a physical 
nanodevice may be synthesized experimentally using a 
combination of the first two prescriptions. 

The structure of the nanodevice in Fig. [l] in partic¬ 
ular seems like this type of quantum dragon should be 
amenable to experimental synthesis. The metal polonium 
(Po) has a simple cubic lattice structure, and hence Fig. [1] 
can be viewed as a nano-crystal of Po. Furthermore, the 
proofs show that any nano-crystal of Po can be connected 
to be a quantum dragon, particularly with the end slices 
as (100) faces as in Fig. [TJ Of course, given the half-life 
of 208 Po is about 2.9 y and that of 209 Po is about 125 y, 
a perfect nano-crystal of Po will only survive so long be¬ 
fore nuclear decay creates defects. It should be possible 
to connect homogeneous leads to a Po nano-crystal, and 
the search for electrical conductivity showing dragon seg¬ 
ments may be helped by shaped electric potentials in the 
nanocrystal. 

The proofs of existence of quantum dragon segements 
here are only for leads of a single channel, and for ho¬ 
mogeneous leads. Furthermore, the proofs are for the 
single-band tight binding model. It may be possible to 
extend these proofs to multi-channel leads, to other more 
complicated leads, and to more realistic band models. 
Of particular interest would be to try to find quantum 
dragons in face-centered cubic single crystal nanodevices. 
This paper has shown how ubiquitous quantum dragons 
are for the simplest cases, but lends hope to their ex¬ 
istence in more complicated nanodevices connected to 
more complicated leads. 

It is anticipated that quantum dragons will at least 
have similar technological applications as do ballistic 
electron propagation devices [g, Effij - EIj . Previously 
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all known electron propagation with complete electron 
transmission, T(E)= 1, were for homogeneous nanode¬ 
vices. Whether quantum dragons that are inhomoge¬ 
neous will enable additional technological applications is 
an active topic of research. 
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APPENDIX 

Herein details describing the four figures are given. Al¬ 
though the figures are somewhat schematic, being only 
four examples of quantum dragons, it is informative nev¬ 
ertheless to give details of their construction. The length 
units are relative since the figures are schematic, but due 
to the atomic nature of the nanodevice, are expected to 
be about a nanometer. All figures were made in Mathe- 
matica [29j. 

Figure 1. Fig. [l] shows the construction of a quan¬ 
tum dragon nanodevice using prescription 1 for a nano¬ 
crystal cut from a simple cubic lattice. The slice in 
Fig. 1(a) is made by keeping only sites from a square 
lattice within an ellipse (green). The square lattice has 
lattice spacing unity, and each atom is a (red) sphere 
of radius 0.5. Both nearest-neighbor (which cannot be 
seen because the spheres touch) and next-nearest hop¬ 
ping terms (blue cylinders) are present. Only the ra=76 
atoms within an ellipse centered at the randomly chosen 
point (0.539, 0.253) with axis along x of length 12 and 
along y of length 2 are used in the slice. The number of 
non-zero hopping terms (number of bonds) in the slice 
is 236. Fig. 1(b) shows the connections between the lead 
(white sphere) and the first slice of the nanodevice, with 
the slice being identical to that of Fig. 1(a). The radius 
of the lead-device (orange) cylinders are proportional to 
the strength of a particular lead-device hopping strength 
s w . The lead-device hopping values are found by find¬ 
ing a vector w\ which is an eigenvector of A with all 
non-positive elements, which is guaranteed to exist by 
the Perron-Frobenius theorem. The input lead atom is 
placed at the CM given by the Wj , but positioned above 
the first slice. Fig. 1(c) shows the completed nanodevice. 
Only two atoms (white spheres) in each semi-infinite lead 
are shown, and for clarity the lead bonds (cyan cylin¬ 
ders) are plotted at a distance of 1.5. The end atom of 
each semi-infinite lead is plotted to be separated from 
the blob by 1.5. There are £=20 slices plotted, separated 
by a distance of 1.1. In Fig. [TJ since ra=76, quantum 
dragon segments exist on a 154-dimensional ‘surface’ in 
3079-dimensional space of all tight-binding parameters. 


Figure 2. Fig. [2] shows the construction of a quan¬ 
tum dragon nanodevice using prescription 1 for a ran¬ 
domly constructed slice, (a) The arrangement of the 
atoms (red spheres) in the slice is constructed by choos¬ 
ing an ellipse (green curve) centered at the origin with an 
x-axis equal to 12 and a y -axis equal to 2. Fifty atoms 
are placed with uniform probability within the ellipse. 
The atoms are assumed to have a hard radius equal to 
0.5, and hence cannot be placed closer than a center-to- 
center distance of djj/= 1.0. The intra-slice bonds are 
shown by (blue) cylinders, with bonds placed between 
any atoms with djj/< 2.5. The strength of the intra-slice 
bonds are given by a linear relationship in the center- 
to-center distance, with bond strength unity for djj >=1 
and zero for djj'= 2.5. (b) The connection of a slice to 
a lead atom (white sphere) in order to satisfy prescrip¬ 
tion 1 is shown. The lead-device hopping strengths are 
calculated by finding a vector w\ which is an eigenvec¬ 
tor of the intra-slice A with all non-positive elements. 
The vector w\ is unique and guaranteed to exist by the 
Perron-Frobenius theorem. The strength of the lead-slice 
bonds are proportional to the radii of the (orange) cylin¬ 
ders. The input lead atom is positioned at the CM given 
by the hopping parameters Wj, but above the plane of 
the first slice, (c) The complete nanodevice, here com¬ 
posed of £=8 identical slices. The transmission can be 
calculated either from Eq. m which requires finding 
the inverse of a (£ra-f2) x (£m+2) =402x402 matrix or 
from Eq. dH which requires finding the inverse of a 
(£+2)x(£+2)=10xl0 matrix. Fig. [2] has m=50, dragon 
segments exist on a 102-dimensional ‘surface’ in the 1376- 
dimensional space of all tight-binding parameters. 

Figure 3. Fig. [3] shows the construction of a slice of 
a quantum dragon nanodevice using prescription 2. The 
lattice spacing is set to one for both the green and blue 
square lattices, (a) The (green sphere) atoms are within 
a generalized ellipse with the equation 



with the center randomly set. All 34 (green) atoms are 
inside this (green) generalized ellipse but not inside the 
blue circle of radius 3 centered at the point (2,2). The 
26 (blue) atoms are inside this (blue) circle. The square 
lattice of these (blue) atoms is randomly offset from the 
center. Any of the (blue) atoms inside the circle that 
would be at a distance less than unity to a (green) atom 
are not included. The slice has m=60 atoms. Green- 
to-green atom bonds (50 bonds) and blue-to-blue atom 
bonds (40 bonds) are only between atoms at a distance 
of unity (nearest neighbors). Green-to-blue atoms bonds 
are inversely proportional to their length, and are be¬ 
tween the atoms at a distance of less than two, giving 26 
such bonds. Therefore there is a total of 116 intra-slice 
bonds (orange cylinders), (b) The yellow sphere shows 
the position of the lead atom in order for it to be at the 
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center of mass [located at (2.08, -0.33)] of the slice for the 
eigenvector of A (as in prescription 1). Instead, the lead 
atom is chosen to be at the point (1,1,2), given by the 
(black) sphere, while the slice is in the plane z= 0. Every 
atom in the slice is connected to the lead atom with a 
strength (black cylinders) chosen to be proportional to 
r~ 6 , with rj the distance between the lead atom and the 
j th atom in the slice, (c) The required electric poten¬ 
tial Vj at every atom site j in the slice is shown (cyan 
cuboid) in order to make the slice satisfy the mapping 
of Eq. (|SJ). For visual reasons, these potentials are all 
shifted by the same amount in order to make them all 
non-negative. For a quantum dragon, typically the re¬ 
quired Vj may be of different signs for different atoms. 
In Fig.[3j since m=60, dragons live on a 122-dimensional 
‘surface’ in the tight-binding 1951-dimensional space. 

Figure 4. Fig.0]shows the construction of a quantum 
dragon nanodevice using prescription 3 for completely 
random inhomogeneous slices. The figure has £=4 slices. 
The ellipse for the random placement of atoms in each 
slice has the x ax i S in Table 1 and y ayi j s =2. The spheres 
have a hard core radius t>, and the intra-slice bonds are 
for any atom pairs in the same slice with a distance less 
than 4. Table 1 shows the specific values for each slice. 
The fixed intra-slice bond strengths are proportional to 
the width of the intra-slice (blue) cylinders, and were cho¬ 
sen as a linear function of distance, with width one for 
djj '=1 and width zero for djj'= 4. The diameters of the 
cylinders of the inter-slice bonds, and of the lead-device 
bonds, are proportional to the tuned bond strengths re¬ 
quired for the device to be a quantum dragon. The input 
(output) lead atoms are placed at the CM of the hopping 
parameters (wij) that connect the input (output) 
lead to the atoms in the first (last or ^ th ) slice, but just 
below (above) the first (last) slice. In Fig. [4j since £=4 
and the rrik are listed in Table 1, the quantum dragons ex¬ 
ist on a 270-dimensional ‘surface’ in the 652-dimensional 
space of all tight-binding parameters. 


TABLE I. Parameters related to Fig. [4] 


Slice # 
k 

m k 

axis 

# intra-slice 
bonds 

1 

8 

2.0 

28 

2 

12 

8.0 

27 

3 

10 

4.0 

33 

4 

14 

10.0 

38 


REFERENCES 


[1] M.A. Novotny, “Energy-independent total quantum 
transmission of electrons through nanodevices with cor¬ 
related disorder”, Phys. Rev. B 90, 165103 [14 pages] 
(2014) 

[2] R. Landauer, “Spatial variation of currents and fields due 
to localized scatters in metallic conduction”, IBM J. Re¬ 
search and Development 1, 223-231 (1957). 

[3] S. Datta, Electronic Transport in Mesoscopic Systems 
(Cambridge University Press, Cambridge, UK, 1995). 

[4] D.K. Ferry and S.M. Goodnick, Transport in Nanos¬ 
tructures (Cambridge University Press, Cambridge, UK, 
1997). 

[5] S. Datta, Quantum Transport: Atom to Transistor 
(Cambridge University Press, Cambridge, UK, 2005). 

[6] G.W. Hanson, Fundamentals of Nanoelectronics , 
(Prentice-Hall, Englewood Cliffs, NJ, 2008). 

[7] Y.V. Nazarov and Y.M. Blanter, Quantum Transport 
(Cambridge University Press, Cambridge, UK, 2009). 

[8] H.-S. Wong and D. Akinwande, Carbon Nanotube and 
Graphene Device Physics (Cambridge University Press, 
Cambridge, UK, 2011). 

[9] Shankar R., Principles of Quantum Mechanics , Second 
Edition (Plenum Press, London, 1994). 

[10] J.W. Mintmire, B.I. Dunlap, and C.T. White, “Are 
fullerene tubules metallic?”, Phys. Rev. Lett. 68, 631- 
634 (1992). 

[11] N. Hamada, S.I. Sawada, and A. Oshiyama, “New one¬ 
dimensional conductors: Graphitic microtubules”, Phys. 
Rev. Lett. 68, 1579-1581 (1992). 

[12] S. Frank, P. Poncharal, Z.L. Wang, and W.A. de Heer, 
“Carbon nanotube quantum resistors”, Science 280, 
1744-1746 (1998). 

[13] M.S. Dresselhaus, G. Dresselhaus, J.C. Charlier, and E. 
Hernandez, “Electronic, thermal and mechanical proper¬ 
ties of carbon nanotubes”, Phil. Trans. R. Soc. Lond. A 
362, 2065-2098 (2004). 

[14] S.J. Tans, M.H. Devoret, H. Dai, A. Thess, R.E. Smalley, 
L.J. Geerligs, and C. Dekker, “Individual single-wall car- 
gon nanotubes as quantum wires”, Nature 386, 474-477 
(1997). 

[15] C. Berger, Y. Yi, Z.L. Wang, and W.A. de Heer, “Multi- 
walled carbon nanotubes are ballistic conductors at room 
temperature”, Appl. Phys. A 74, 363-365 (2002). 

[16] D. Daboul, I. Chang, and A. Aharony, “Series expan¬ 
sion study of quantum percolation on the square lattice”, 
Euro. J. Phys. B 16, 303-316 (2000). 

[17] M.F. Islam and H. Nakanishi, “Localization- 
delocalization transition in a two-dimensional quantum 
percolation”, Phys. Rev. E 77, 061109 [9 pages] (2008). 

[18] E. Cuansing and J.-S. Wang, “Quantum transport in 
honeycomb lattice ribbons with armchair and zigzag 
edges coupled to semi-infinite linear chain leads”, Euro. 
Phys. J. B 69, 505-513 (2009). 

[19] S. Boettcher, C. Varghese, and M.A. Novotny, “Quantum 
transport through hierarchical structures”, Phys. Rev. E 
83, 041106 [12 pages] (2011). 










14 


[20] Z. Lin and Y. Liu, “Electronic transport properties of the 
Bethe lattices”, Phys. Lett. A 320, 70-80 (2003). 

[21] C. Varghese and M.A. Novotny, “Quantum transport 
through fully connected Bethe lattices”, Inter. J. Mod. 
Phys. C 23, 1240010 [10 pages] (2012). 

[22] R.G. Busacker and T.L. Saaty, Finite Graphs and Net¬ 
works: An Introduction with Application (McGraw-Hill, 
New York, 1965). 

[23] G. Chartrand and P. Zhang, A First Course in Graph 
Theory (Dover Publications, New York, 2012). 

[24] M. Marcus and H. Mine, Matrix Theory and Matrix In¬ 
equalities (Dover Scientific, New York, 1992) from the 
1964 original edition. 

[25] A. Berman and R.J. Plemmons, Nonnegative Matrices in 
the Mathematical Sciences (Academic Press, New York, 
1979). 


[26] A. Javey, J. Guo, Q. Wang, M. Lundstrom, and H. Dai, 
“Ballistic carbon nanotube field effect transistors”, Na¬ 
ture 424, 654-657 (2003). 

[27] J. Wu, L. Xie, G. Hong, H. En Lim, B. Thendie, Y. Miy- 
ata, H. Shinohara, and H. Dai, “Short channel field-effect 
transistors from highly enriched semiconducting carbon 
nanotubes”, Nano Research 5, 388-394 (2012). 

[28] S.H. Kim, W. Song, M.W. Jung, M-A. Kang, K. Kim, S- 
J. Chang, S.S. Lee, J. Lim, J. Hwang, S. Myung, and K-S. 
An, “Carbon nanotube and graphene hybrid thin film for 
transparent electrodes and field effect transistors”, Adv. 
Mater. 26, 4247-4252 (2014). 

[29] Wolfram Research, Inc., Mathematica, Version 9.0, 
Champaign, IL (2012). 



